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Abstract We re-analyze Nuclear Magnetic Resonance (NMR) spectra observed at low tem- 
peratures and high magnetic fields in the field-induced B -phase of CeCoIns. The NMR 
spectra are consistent with incommensurate antiferromagnetic order of the Ce magnetic mo- 
ments. However, we find that the spectra of the In(2) sites depend critically on the direction 
of the ordered moments, the ordering wavevector and the symmetry of the hyperfine cou- 
pling to the Ce spins. Assuming isotropic hyperfine coupling, the NMR spectra observed 
for H || [100] are consistent with magnetic order with wavevector Q = 7z(^^-,-,-) and 
Ce moments ordered antiferromagnetically along the [100] direction in real space. If the 
hyperfine coupling has dipolar symmetry, then the NMR spectra require Ce moments along 
the [001] direction. The dipolar scenario is also consistent with recent neutron scattering 
measurements that find an ordered moment of 0.15^ along [001] and Q n = it(^- 1 -) 
with incommensuration 8 =0.12 for field H | [110]. Using these parameters, we find that a 
hyperfine field with dipolar contribution is consistent with findings from both experiments. 
We speculate that the B phase of CeCoIns represents an intrinsic phase of modulated super- 
conductivity and antiferromagnetism that can only emerge in a highly clean system. 
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1 Introduction 

The heavy-fermion superconductor CeCoIns has attracted considerable attention since its 
discovery in 2001. 1 Not only does this unconventional d-wave superconductor exhibit non- 
Fermi liquid behavior associated with proximity to a proposed quantum critical point, but it 
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Fig. 1 (Color online) The phase diagram of CeCoIns 
in high field as determined by specific heat.— Solid 
points represent second order phase transitions and 
open points are first order transitions. The solid blue 
squares are the points at which the spectra in Fig. [3] 
were obtained. 




Fig. 2 (Color online) The unit cell of CeCoIns . The 
Ce atoms (yellow) sit at the eight corners. The In(l) 
atoms sit in the center of the top and bottom faces (or- 
ange). The Co atoms are grey and the In(2) atoms are 
green. For the field oriented in the ab plane, there are 
two inequivalent In(2) atoms, depending on whether 
the field is parallel (In(2a) or perpendicular (In(2b)) 
to the unit cell face. 



is unique among the heavy-fermion superconductors in that it also exhibits a new thermo- 
dynamic phase (B phase) that exists only within the superconducting phase near H c 2-^ Ini- 
tially this B phase was identified as the elusive Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) 
superconducting phase first predicted to exist in Pauli-limited superconductors over 40 years 
ago. 4 ] 5 ] 6 ] 7 ! 8 In fact, recent NMR work by Young and coworkers 9 identified the presence 
of incommensurate antiferromagnetic order in the B phase in contrast to the standard pre- 
dictions for the FFLO phase. 10,11 Signatures of magnetism were also seen in other NMR 
experiments. 12 i 13 The NMR spectra of the In sites in the B phase do not reveal additional 
paramagnetic resonances that might be associated with macroscopic phase separation, but 
rather are consistent with superconducting order coexisting with a homogeneous modulation 
of antiferromagnetic order. 9 Despite initial arguments to the contrary,— recent neutron scat- 
tering results by Kenzelmann and coworkers now provide conclusive proof for long-range 
static incommensurate antiferromagnetic order.— 

The antiferromagnetism in CeCoIng was first identified by Young et al. due to the pres- 
ence of a broad spectrum observed at the In(2) sites in this material (see Figs. [2] and [3). 
The In(l) and Co sites, in contrast, showed no splitting. Young et al. pointed out that these 
observations place constraints on the possible magnetic structure, but do not uniquely iden- 
tify the structure. They proposed a minimal model where the magnetic structure consists of 
ordered local Ce spins with moments So along the applied magnetic field direction (along 
[100]), with an ordering wavevector of the form Q = n{ iii , ~, i). The structure of the NMR 
spectra revealed the incommensurate nature, but the value of the modulation 8 remained un- 
determined since the hyperfine field at the In(2) site depends on the product of the size of the 
ordered moment and the incommensuration. In the neutron diffraction experiment, Kenzel- 
mann et al. oriented the field along the [1 10] direction, and observed Q„ = 7r(i±^, ijp., |) 
with 8 = 0.12 and moments along [001]. A crucial observation was that 8 is independent 
of the applied field in the B phase, in contrast to the predictions for the FFLO phase. By 
proposing a Ginzburg-Landau model for the coupling of antiferromagnetism and supercon- 
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Fig. 3 (Color online) Fixed field NMR spectra at 11.1 T in CeCoInj showing how the In(l), In(2) and Co 
sites evolve from the normal state (top) to the B phase (bottom), adapted from Young et al— 



ductivity, they showed that the superconducting order parameter in the B phase acquires a 
component with finite momentum because of the strong coupling between the incommen- 
surate magnetism and superconductivity. The neutron data confirm the NMR observation 
that this exotic state disappears immediately above H c 2, where the system returns to a fully 
homogeneous normal phase, yet with strong deviations from conventional Fermi liquid the- 
ory.— 

A priori, the NMR and neutron scattering results suggest different magnetic structures. 
In order to address this discrepancy, we investigate several possible magnetic structures al- 
lowed by the NMR results. The neutron diffraction results suggest that the applied field, H, 
the moments, So, and the incommensuration wavevector Q, = - (8,8,0) are all mutually 
orthogonal. In contrast, the proposed NMR scenario suggested So || H || Qo. As we show 
below, this scenario is the most likely for isotropic transferred hyperfine couplings between 
the Ce spins and the In(2) nuclei. On the other hand, if the coupling tensors are anisotropic, 
then other magnetic structures are possible as argued by Koutroulakis et al. 16 i 17 If the hy- 
perfine tensor has purely dipolar symmetry, then we find that for H | [100] the most likely 
magnetic structures satisfy So 1 1 [001] and Q, ± So- The spectra of the In(l) and In(2) differ 
slightly depending on the orientation of Q; in the plane, but the data are most consistent 
with Qj 1 1 [010] or Q,- 1 1 [110]. The anisotropic coupling scenario offers a picture that is both 
physically more reasonable and consistent with neutron diffraction observations of Q, _L So 
and H 1 S , though for fields H 1 1 [110]. 
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2 Analysis 

2.1 Spectra 

In order to explain the broad double-peak structure of the In(2) spectrum in Fig. [3] there must 
be a distribution of local fields both parallel and antiparallel to the applied field Ho||[100] 
with values ranging up to 1.3 kOe.~ The resonance frequency is given by / = y|Ho + Uhf\ + 
/g, where y is the gyromagnetic ratio and /g is the contribution from the quadrupolar in- 
teraction at the nucleus. Since the electric field gradient (EFG) at all three nuclear sites is 
unaffected by the onset of superconductivity or magnetism, the dramatic line broadening 
effects observed in the B phase can be attributed entirely to the onset of the static hyperfine 
field, H/,y. /g is a temperature independent constant that depends on the particular site and 
we will not address it further. Experimentally, we find no significant broadening in the B 
phase for the In(l) or the Co, but a broad, double-peak spectrum for the In(2a) site (previ- 
ously referred to as In(2)n, the In(2) site on the unit cell face that lies parallel to the field). 
Independent measurements show no significant broadening at the In(2b) (previously referred 
to as In(2)x). These results put stringent constraints on any candidate magnetic structure. 

The double-peak structure observed for the In(2a) arises because there is a distribution 
of local hyperfine fields that lie either parallel or antiparallel to the applied field. If H/,y is 
parallel to Ho and is modulated along a direction f, then f(r) = fo + yh^j- cos(qr), where 
n^f is the magnitude of the modulation, r is the distance along the modulation, and q is the 
wavevector. In this case, the spectrum will then be given by 



\(f)cc\dfldr\- 1 = , * (1) 



f(h° hf ) 2 -(f-f ) 



On the other hand, if Hfty(r) _L Ho then /(r) = W/q + y 2 {h a h y) 2 cos 2 (qr) and the spectrum 
is given by 

&x if) « ; ^-L ■ (2) 

\ff^7^f(K f ) 2 -f 1 -fZ 

These spectra are shown in Fig. [4] and clearly show that for the parallel case, there are two 
double peaks at frequencies both below and above /o, whereas for the perpendicular case, 
there is only a single peak at higher frequency. A key result of the re-analysis of the NMR 
spectra is, in agreement with our earlier analysis, that the In(2a) sites require a hyperfine field 
parallel to the applied magnetic field, H;,/ 1 1 Ho, to account for the broadening and double- 
peak spectrum. On the other side, the In(2b) sites show little or no broadening, consistent 
with a vanishing or perpendicular hyperfine field, i.e., at the In(2b) sites either H/,y = or 
H/,/ -L Hq. 



2.2 Hyperfine couplings 

The hyperfine interaction is given by the Hamiltonian 

^ / = l-A-S(r = 0)+£l-B,-S(r / ), (3) 

i 

where the hyperfine coupling tensor A represents the on-site coupling to an electron spin, S, 
at the nuclear site r = 0, B represents a transferred hyperfine coupling to an electron spin on a 
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Fig. 4 (Color online) The calculated spectra of &\\{f) (blue) and ,3?±(f) (red) assuming /o= 100 MHz and 
yh®f = 1 MHz. For the perpendicular case, the spectrum is only weakly affected by the hyperfine fields, 
whereas for the parallel case it broadens dramatically resulting in a double-peak structure. 



distant (ligand) site at r,. 18 In CeCoIng, these sites are the nearest neighbor Ce 4/ electrons, 
and the sum is over the nearest neighbors. For static ordering of the Ce spins, Eq. (O can be 
re-written as: Jf?/, t = yhl ■ H/,y, where the magnitude and direction of the hyperfine field H/, / 
depend critically on the hyperfine tensors for the particular site and the magnetic structure. 
The tensorial A term represents hyperfine coupling to the itinerant conduction electrons, 
which we will ignore, since we are only concerned with static contributions to H/,/ from the 
static local Ce ordering. 

The transferred hyperfine tensor B is generally not diagonal in the crystal axis basis,— 
and may be written as the sum of isotropic and dipolar contributions. 20 To lowest order, the 
tensor can be approximated by a scalar (isotropic) interaction, since the transferred hyperfine 
interaction is typically at least one order of magnitude greater than the direct dipolar inter- 
action. However, in the CeCoIng compound, there is evidence that the hyperfine interaction 
is not purely isotropic, and therefore we must consider dipolar symmetries as well. Indeed, 
the magnitude of the dipolar portion is found to be enhanced by the delocalized nature of 
the electrons in the solid. 21 Therefore, we write the hyperfine fields at the ligand sites as: 

4 

Hv(r) = £B r S(r + r,)/yfi atln(l) (4) 

i=i 

2 

H ft/ (r) = £B r S(r + rj)/yfi atCo (5) 

;'=i 

2 

H v (r) = £B r S(r+rk)/7fi atln(2a) (6) 

jt=i 

2 

H v (r) = £B i -S(r+n)/7», atln(2b) (7) 

z=i 

where r; = (±f ,±|,0) for the In(l) nearest neighbor Ce sites, rj = (0,0, ±|) for the Co 
nearest neighbor Ce sites, and r/ = (±|,0,zo) for the In(2a) nearest neighbor Ce sites, and 
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tk = (0,±§,zo) for the In(2b) nearest neighbor Ce sites. The couplings 
where: 

/l 00\ 



10 (8) 
00 1/ 



and 

B I 2x 2 - y 2 - z 2 3xy 3xz \ 

B d ip = ^M 3xy -x 2 + 2y 2 -z 2 3yz (9) 
r \ 3xz 3yz -x 2 -y 2 + 2z 2 J 

are evaluated for each site. Here r = (x,y,z) is the vector joining the particular site to the Ce 
atom in question. 



Table 1 The hyperfine fields at the In(l) and In(2) sites for various magnetic structures. Here we assumed 
that the applied magnetic field is Ho || [100], the moment is 0.15 iig and its modulation is S = 0.12. The 
column at the far right indicates magnetic structures in agreement with the NMR spectra. 
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dip 
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no 


(1.3) 


[100] 
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Hhf 
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yes (e) 



2.3 Magnetic structure 

The magnetic structure is given by S = Socos[(Qo + Q/) -r], where the antiferromagnetic 
wavevector Qo = (f , f , f) is commensurate with the lattice. Neutron diffraction reports an 
incommensurate wavevector Q; = -(8,8,0) with spatial modulation \/2a/8 « 12a w 5.4 
nm in the ab plane, and ordered moment So at the Ce site. This modulation is significantly 
shorter than the average inter-vortex distance of ~ 14 nm in a field of 10 T. 13 So it does not 
support a picture of overlapping extended states from the vortex cores leading to magnetic 
ordering. Similarly this specific Q, would be inconsistent with a field along [100], if it was 
related to the alignment of vortices along the [100] direction. In our re-analysis of the NMR 
spectra, we consider four cases: (1) Q,- || [100], (2) Q, || [010], and (3) Q, || [110], (4) 
Q,- || [110]. Case (1) was proposed by Young et al.— for NMR measurements under the 
condition Hq 1 1 [100]. Case (2) should be equally likely as case (1) because of the tetragonal 



Fig. 5 (Color online) Magnetic structure and hyper- 
fine fields for isotropic hyperfine couplings to the 
In(2) (case (1.1)). The in-plane tetragonal structure 
is outlined in gray. The Ce atoms are yellow, and 
their moments are indicated by red arrows pointing 
along [100]. The In(l) atoms are orange, and the Co 
are not shown. The In(2a) are green and the In(2b) 
are blue. The hyperfine fields at the In(2a) sites are 
indicated by blue arrows. Here the hyperfine fields 
vanish at the In(l), Co and the In(2b) sites; the di- 
rection of Ho is shown by the black arrow. 



Fig. 6 (Color online) Magnetic structure and hy- 
perfine fields for dipolar hyperfine couplings to the 
In(2) (case (2.4)). Same notation as in Fig. \5\ The 
Ce moments point along [001]. Here the hyperfine 
fields vanish at the In( l) and Co sites, but not at the 
In(2) sites; the direction of Ho is shown by the black 
arrow. 



symmetry of the crystal structure. Case (3) was proposed by Kenzelmann et al. — for neutron 
diffraction measurements under the condition Ho || [110], and case (4) should be equally 
likely as case (3). We have calculated the hyperfine fields for each of these cases for both 
purely isotropic and purely dipolar couplings, for Ho 1 1 [100] with moments along both [100] 
and [001], and the results are summarized in Table (1). 

The cases that are most consistent with the NMR observations are (1.1) Q; || [100], 
So || [100] and isotropic coupling, (1.4) Q, || [100], So || [001] and dipolar coupling, and 
(2.4) Q,- 1 1 [010], So| | [001] and dipolar coupling. Cases (3.4) and (4.4) are consistent with the 
In(2) spectra, but give rise to an internal field at the In(l) site that at first sight is inconsistent 
with experiment. We will discuss this in more detail below in the discussion section. Figures 
|5]and|6]show the magnetic structure and hyperfine fields for cases (1.1) and (2.4). Case (1.1) 
is identical to the one originally proposed by Young and coworkers, 9 which most likely will 
not minimize the magnetic contribution to the free energy, as the moments are either parallel 
or antiparallel to the applied field. Cases (1.4), (2.4), (3.4) and (4.4), in which the moments 
are perpendicular to the applied field, are physically more reasonable for antiferromagnetic 
ordering and agree with the neutron diffraction results of Q, ± So and Hq _L So.— 



3 Discussion 

3.1 Hyperfine constants from the Knight shift 

In CeCoIns measurements of the Knight shift in the normal state show that the hyperfine 
couplings for the In(2a) site (r,t = (±|,0,zo)) for the applied field along (100), (010) and 
(001) are 10.3 kOe/,u B , 0.0 kOe//i B , and 32.4 kOe/,u B , respectively. 22 For this site, Equations 
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Fig. 7 (Color online) Real-space map (15 x 15 unit 
cells) of the hyperfine field at the In(2a) (Eq. 6) (up- 
per row) and the In(2b) (Eq. 7) (middle row) in the 
ab plane. Shown along the horizontal are the compo- 
nents of the hyperfine field along the (x, y, z) direc- 
tions in (red, green blue) shading, for case (2.4) with 
Qo || [010] and So || [001]. Black corresponds to zero 
hyperfine field. The lower row shows the spin density 
along the [001] direction (cyan), and the histogram of 
resonant frequencies for the In(2a) (orange) and the 
In(2b) (yellow). 




Fig. 8 (Color online) The same notation as in Fig. [7] 
but for case (3.4) with Qo 1 1 [1 10] and S 1 1 [001]. The 
corresponding histogram of resonant frequencies is in 
better agreement with experiment. 



© and © yield: 

K a = {2B Lm +B dip {\ -3cos2G z ))x a (10) 

K h = (2B iso -2B dip )x b (11) 

K c = (2B iso +B dip (l+3cos28 z ))Xc, (12) 

where % a is the susceptibility in the a direction. Using the experimental numbers, 22 we find 
Biso = B dip =7.1 kOe//i B and Q z = 29°. The difference between this angle and that of the 
crystal structure (6 Z = 45°) probably is related to details of the bonding of the In 4p orbitals, 
and will need further investigations. 



3.2 Hyperfine fields in the B phase 



For cases (1.4), (2.4), (3.4) and (4.4) the hyperfine field at the In(2a) and In(2b) sites varies 
spatially ~ cos[(Qo +Q,) ■ r] with modulus: 

B hf (2a) = 5 5 rfl> [3sin(2a)cos(^)a + (l-cos(2a)sin(^)c] (13) 
K hf (2b) = S o B di p[3sin(20 z )cos(^)b+(l-cos(20 z )sin(^)e] (14) 

In each case, the hyperfine field at the In(2a) site oscillates along the modulation direction 
with a component along [100] and the resulting spectrum is described by Figs . [3] [7] and [8] 
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Using the values 8 X = 8 y « 0.12 and So « 0.15^ as reported by Kenzelmann et al.,— we 
find h®j ks 2.6 kOe, which is about twice the experimental value of 1.3 kOe. The difference 
may be related to uncertainties in the hyperfine coupling itself 18 or changes in the magnetic 
structure for the field along [100]. Note that the modulation and ordered moment may differ 
for field oriented along [100], in which case Eq. d 1 3 b will give a different value. In fact the 
NMR result is consistent with the locus of points given by <5 A = <5 t o cos "'(Sg/So-l), where 
5,0 = 0.3 and Sg = 0.14^. 

For cases (3.4) and (4.4), where the modulation is along [110] or [110], the hyperfine 
field at the In(l) site does not cancel but has a component along the [001] direction. This 
field can give rise to a minor shift and/or broadening of the In(l) line. The spectra (Fig. 
[3} clearly show that the In(l) line shifts and is only slightly broadened. However, the shift 
and broadening may come from the onset of spin shift suppression in the superconducting 
state and the presence of superconducting vortices. Therefore, we cannot distinguish the 
presence of a hyperfine field from the antiferromagnetic structure at the In(l) site within 
experimental error. As seen in Figs. [7] and [8] the calculated spectra for case (3.4) is closer 
to the experimental one. We speculate that the true magnetic structure for the field Ho along 
[100] is best described either by case (2.4), (3.4) or (4.4) with Ho _L So, which will minimize 
the free energy of the antiferromagnet. In each case the hyperfine field at the In(2a) and 
In(2b) have components perpendicular to Ho, but these components only give rise to small 
shifts of the resonant frequency that are difficult to distinguish from the Knight shift. The 
crucial element is that the hyperfine field at the In(2a) is along Ho. 

3.3 Nature of the B phase 

The fact that the antiferromagnetism exists only in field and only within the superconduct- 
ing phase indicates a strong coupling between these order parameters. Kenzelmann and 
coauthors 15 analyzed the symmetry of the superconducting state for such a coupling and 
concluded that the superconducting order parameter, Aq acquires the finite momentum Q 
of the antiferromagnetic order parameter, Mq. This corresponds to a modulation of the or- 
der parameter in real space that presumably is out-of-phase with the antiferromagnetism. In 
other words, the antiferromagnetic order is maximum at the nodes of Aq. Since then vari- 
ous microscopic models have been proposed to explain the field induced antiferromagnetic 
order. 23,24 Curiously, this scenario is similar to that observed in the ferropnictide SrFe2As2 
under pressure. 2 - In this compound, a novel hybrid state of coexisting superconductivity 
and antiferromagnetism emerges above 5 GPa. We speculate that these two novel states 
may in fact be the same. Although superconductivity and antiferromagnetism are known to 
coexist inhomogeneously in a number of doped high T c , heavy fermion, and ferropnictide 
systems, the highly clean undoped CeCoIns and SrFe2As2 materials support the emergence 
of this fragile but intrinsic thermodynamic phase of modulated antiferromagnetism and su- 
perconductivity. In CeCoIns, this phase is quickly destroyed by doping and is replaced by 
a commensurate order at zero field for sufficiently high Cd doping. 26 i 27 i 28 - 29 Clearly many 
questions about this new state of matter remain unexplained, such as the driving mecha- 
nism^), the origin of the incommensurate wavevector, and the nature of the excitations. 

4 Conclusions 

In summary, we have shown that both NMR and neutron diffraction measurements in the 
field-induced B phase of CeCoIns are consistent with magnetic structures where Q, _L So 
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and Ho _!_ So- The incommensurate modulation Q, lies possibly along either [010] or [110] 
direction for magnetic field pointing along [100] in real space. Tetragonal equivalent direc- 
tions for Q,-, [100] and [110], are also possible. Based on our analysis of the NMR spectra, 
we speculate that the B phase of CeCoIns represents an intrinsic phase of modulated super- 
conductivity and antiferromagnetism that can only emerge in a highly clean system. Further 
NMR and neutron diffraction measurements are necessary for the same field orientations to 
unravel the origin of the field-induced antiferromagnetic structure. 
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